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Abstract. Lens time delays are a powerful probe of cosmology, provided that the gravitational 
potential of the main deflector can be modeled with sufficient precision. Recent work has shown 
that this can be achieved by detailed modeling of the host galaxies of lensed quasars, which appear 
as “Einstein Rings” in high resolution images. The distortion of these arcs and counter-arcs, as 
measured over a large number of pixels, provides tight constraints on the difference between the 
gravitational potential between the quasar image positions, and thus on cosmology in combination 
with the measured time delay. We carry out a systematic exploration of the high resolution imaging 
required to exploit the thousands of lensed quasars that will be discovered by current and upcoming 
surveys with the next decade. Specifically, we simulate realistic lens systems as imaged by the Hubble 
Space Telescope (HST), James Webb Space Telescope (JWST), and ground based adaptive optics 
images taken with Keck or the Thirty Meter Telescope (TMT). We compare the performance of 
these pointed observations with that of images taken by the Euclid (VIS), Wide-Eield Infrared Survey 
Telescope (WEIRST) and Large Synoptic Survey Telescope (LSST) surveys. We use as our metric the 
precision with which the slope y' of the total mass density profile ptot ^ for the main deflector 
can be measured. Ideally, we require that the statistical error on y' be less than 0.02, such that it 
is subdominant to other sources of random and systematic uncertainties. We find that survey data 
will likely have sufficient depth and resolution to meet the target only for the brighter gravitational 
lens systems, comparable to those discovered by the SDSS survey. Eor fainter systems, that will be 
discovered by current and future surveys, targeted follow-up will be required. However, the exposure 
time required with upcoming facilitites such as JWST, the Keck Next Generation Adaptive Optics 
System, and TMT, will only be of order a few minutes per system, thus making the follow-up of 
hundreds of systems a practical and efficient cosmological probe. 


Contents 


1 Introduction 1 

2 Summary of simulated instrumental setups 3 

3 The Lens Sample 5 

3.1 Sample Selection 5 

3.2 Parameters of the Sample of Mock Lenses 6 

4 Description of the inference process 7 

4.1 Image generating process 7 

4.1.1 Light Model 7 

4.1.2 Mass Model 8 

4.1.3 Point Source Model 8 

4.1.4 Image Noise 8 

4.2 Inferring the parameters 9 

4.3 Exposure times 9 

5 Results 10 

6 Summary 11 


1 Introduction 

In the past few years, gravitational time delays (Refsdal 1964 [1]) have emerged as a powerful and 
cost effective cosmological probe. Studies based on blind analysis have shown that a single system 
consisting of multiple images of a background quasar (at redshift Zs) lensed by a foreground massive 
elliptical galaxy at redshift Zd can be used to measure the so-called time delay distance with an uncer¬ 
tainty of 5-7% (Suyu et al. 2010 [2], 2013 [3], 2014 [4]). The time delay distance D/^t gives a direct 
measurement of the Hubble constant and allows one to break some of the main degeneracies in the 
interpretation of cosmic microwave background data, thus providing tight constraints on parameters 
such as curvature and dark energy equation of state (Coe & Moustakas 2009 [5], Linder 2011 [6], 
Weinberg 2013 [7], Suyu 2012 [8], Suyu et al. 2014 [4]). The time delay distance measurement of Hq 
is comparable in terms of information content to that obtained via the cosmic distance ladder (Riess 
et al. 2011 [9], Freedman et al. 2012 [10]), in that is based on absolute distances and completely 
independent of the properties of the early universe (z > Zs). Importantly, time delay distances are 
independent of the local distance ladder, and thus provide a crucial test of any potential systematic 
uncertainties. Furthermore, being independent, the cosmic distance ladder and time delay distance 
constraints on Hq can be statistically combined for additional gains in precision. 

From an observational point of view, the attainment of precise (i.e. small random errors) and 
accurate (small systematic errors) time delay measurements relies on a few important ingredients. 
First, monitoring of the lensed quasars is required to obtain time delays with a few percent uncer¬ 
tainties. Dedicated monitoring campaigns have shown that this precision is achievable with l-2m 
class telescopes at present time (e.g. from COMOGRAIL; Tewes et al. 2013 [11]), or in the radio 
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(Fassnacht et al. 2002 [12]). In the future, the exploitation of larger samples of lenses (Oguri & Mar¬ 
shall 2010 [13]) will require monitoring campaigns on 2-4m telescopes (Treu et al. 2013 [14]), or 
the deployment of long term high-cadence monitoring efforts, like the LSST (Liao et al. 2015 [15]). 
Second, the spectroscopic redshift of the source and deflector must be measured. For current samples, 
this is typically a relatively straightforward step, requiring short exposures on 4-10m class telescopes 
(Fassnacht & Cohen 1998 [16], Eigenbrod et al. 2007 [17]). Third, the gravitational potential of the 
main deflector needs to be constrained by the data so that the uncertainty on its difference between the 
location of the images is also of order 3%. This goal can be achieved by studying the extended struc¬ 
ture of the lensed quasar host galaxy, and the stellar kinematics of the deflector galaxy (Kochanek 
et al. 2001 [18], Treu & Koopmans 2002 [19], Koopmans et al. 2003 [20], Suyu et al. 2010 [2], 
Suyu et al. 2013 [3], Suyu et al. 2014 [4]). Fourth, the combined effects of the inhomogenous mass 
distribution along the line of sight need to be taken into account. Recent work has shown that the line 
of sight effects can be sufficiently characterized by measuring the properties of galaxies and weak 
lensing signal in the field of the main deflector, and comparing with simulated lines of sight (Suyu et 
al. 2010 [2], 2013 [3], 2014 [4], Greene et al. 2013 [21], Collett et al. 2013 [22]). Fifth, the mass 
sheet degeneracy (Falco et al. 1985 [23]) and its generalizations (Schneider & Sluse 2013, 2014, Xu 
et al. 2015; [24-26]) must be broken either by the incorporation of non-lensing data such as stellar 
velocity dispersion (Treu & Koopmans 2002, Koopmans et al. 2003, Suyu et al. 2014, [4, 19, 20]), or 
by appropriate asssumptions about the mass distribution in the main deflector (Xu et al. 2015 [26]). 

However, whereas current samples have been limited to a few lenses, current and ongoing 
surveys—such as the Dark Energy Survey (DBS), PanSTARRSl, HyperSuprimeCam (HSC)— should 
together discover hundreds of lensed quasars suitable for monitoring and follow-up (Oguri & Mar¬ 
shall 2010 [13]). During the next decade, the Large Synoptic Survey Telescope (LSST), Euclid, and 
the Wide-Field Infrared Survey Telescope (WFIRST) will come online, and enable the detection and 
measurement of an order of magnitude more, with LSST likely to provide high accuracy time delays 
for several hundred systems (Liao et al. 2015 [15]). These large samples will be extremely pow¬ 
erful cosmographic probes (Coe & Moustakas 2009 [5], Linder 2011 [6]), provided that sufficient 
follow-up data are available (Treu et al. 2013 [14]). In this paper we carry out a feasibility study 
for cosmology with future time delay lenses, focusing on the high-resolution imaging follow-up re¬ 
quirements. The multi-pronged nature of time delay lens follow-up makes it natural to follow the 
approach of addressing each component of the follow-up independently. The monitoring and spec¬ 
troscopic follow-up requirements are described, for example, by Eigenbrod et al. (2005; [27]) and 
Linder (2015; [28]). Specifically, we aim to answer the following questions: 

• Will sufficient information be available from the survey data itself, or will dedicated follow-up 
be necessary? 

• If dedicated follow-up is necessary, approximately how much exposure time will be required 
per system, to enable high precision cosmography? 

In order to make the problem tractable and the results general, we need a single simple metric to 
evaluate the quality of imaging data. In short, we need to quantify our ability to measure how the 
deflection angle (i.e. the derivative of the lensing potential) varies between the images, and thus how 
the extended images of the quasar host galaxy are stretched across the image. In practice, our ability 
to constrain the differential magnification will depend on the resolution of the data as well as on the 
number of pixels where the source is detected above a certain signal to noise ratio (and thus on the 
source magnification for fixed instrumental configuration and source intrinsic luminosity). 

We choose to adopt as our metric the slope y' of the total mass density profile of the form ptot 
r~^' as measured by fitting elliptical power law models to the data. This profile is the simplest one that 
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provides a realistic description of galaxy scale lenses (Treu 2010 [29]), and the uncertainty 7 ' has been 
shown to be approximately proportional to the uncertainty on the gravitational potential differences 
and thus the time delay distance, and the Hubble constant (Kochanek 2002 [30], Wucknitz 2002 [31], 
Suyu 2012 [ 8 ]). Intuitively, 7 ' is directly related to the variations of the lensing potential. In fact, if 
y - 2 (the so-called isothermal profile) the deflection angle is constant across the image and therefore 
all the images will appear to have the same radial magnification. If the profile is steeper/shallower, 
radial magnification will vary across the image, thus giving rise to images of different widths. 

Naturally, for a real measurement it is important to explore different profiles and the role of 
the choice of the profile in the uncertainties (see e.g. Suyu et al 2014 and references therein [4]). 
However, our goal is to estimate minimal requirements on the data quality. If the data quality is 
insufficient to constrain 7 ' it will also be insufficient to constrain more flexible models. Likewise, this 
paper is only concerned with precision, i.e. random errors. Systematic errors not considered in this 
paper, such as those arising from incomplete knowledge of the PSF will only increase the error budget 
(e.g., Agnello et al. 2015; Rusu et al. 2015; [32, 33]). Thus, in order to leave room for additional 
sources of uncertainty in the error budget, we set a rather stringent requirement of 0.02 uncertainty 
of y (corresponding approximately to 2 % per system on time delay distance), in tests where the 
mock data are generated using the same model as is used in the inference of the lens parameters. In 
other words, our target corresponds to the requirement that the statistical error arising from the image 
quality be subdominant with respect to those arising from modeling uncertainties, time delays, and 
line of sight effects. 

This paper is organized as follows. In Section 2, we summarize the characteristics of the tele¬ 
scopes and instruments simulated as part of this work. Then, in Section 3, we describe the properties 
of the simulated lenses. Next, in Section 4 we describe the procedure used to carry out the inference. 
Our results are presented in Section 5. Finally, we discuss and summarize our work in Section 6 . 
Throughout this paper, all magnitudes are given in the AB system. Even though our findings are in¬ 
dependent of any cosmology, we adopt a spatially flat ACDM cosmology with Qm = 0.3, = 0.7, 

and the Hubble constant Hq = 70 km s“^ Mpc“^ when calculating distances. 

2 Summary of simulated instrumental setups 

We aim to carry out a systematic exploration of imaging requirements given current and future facil¬ 
ities. For a given lens system, the instrumental setup drives our ability to constrain 7 ' in a number of 
ways. First, the signal-to-noise ratio depends on the exposure-time adopted, sky or background noise 
in a chosen band, and the instrumental readout noise. Second, the pixel-size and the point-spread 
function (PSF) properties determine how finely one can map the system being observed, hence how 
robustly the deflections on either side of the lens can be quantified. 

Here we summarize the main properties of instruments that will be considered in our simula¬ 
tions. In an effort to balance completeness with feasibility, we consider a suite of instruments that 
include ground and space based telescopes, and cadenced surveys where the total exposure cannot be 
set arbitrarily. Some of the choices are representative of a class of telescopes/instruments, for exam¬ 
ple we expect the performance of Keck Adaptive Optics to be a useful guidance for adaptive optics 
systems on other 8 -10m class telescopes, and that of the Thirty Meter Telescope to be a guidance 
for other planned 20-30m class telescopes. For survey telescopes like WFIRST, LSST, and Euclid, 
we only consider the wide field surveys, as the deep surveys cover too small a solid angle to include 
randomly a useful number of lensed quasars. 

Table 1 and figure 1 display the main properties of the telescopes that are needed when generat¬ 
ing images of mock gravitational lenses, taken from the instrument websites. Figure 2 shows the typ- 
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ical PSF of each instrument. In order to facilitate comparisons between instruments and telescopes, 
whenever possible we selected the filter/band for each configuration that is closer to the /-band in the 
optical and the /f'-band in the near infrared ^ A brief description of each setup follows. 

1. We take the Advanced Camera for Surveys (ACS) as our reference optical imager on board 
the Hubble Space Telescope (HST). Its properties are taken from the Space Telescope Sci¬ 
ence Institute HST Exposure Time Calculator.^ In the chosen F814W filter its performance is 
comparable to that of WFC3 (slightly higher sensitivity and coarser pixel size). For simplic¬ 
ity, we neglect charge transfer inefficiency effects, assuming that they are negligible or can be 
corrected to the desired level. The PSF is simulated using the Tiny Tim software (Krist 1993 
[36]). 

2. The Near Infrared Camera (NIRCAM) on board the James Webb Space Telescope (JWST) is 
chosen as the next space based imaging capability. We select the broad F200W filter, where 
the image quality is virtually diffraction limited. We use the instrument properties as given by 
the JWST Exposure Time Calculator.^ 

3. The current and planned adaptive optics systems (hereafter LGSAO & NGAO;[37], [38], re¬ 
spectively) at the W.M.Keck Observatory are chosen to represent current and upcoming AO 
performance on 8-10m class telescopes. We consider the current instrument NIRC2, although 
further gains might be possible with an instrument upgrade in conjuction with the AO system 
upgrade. We also assume conservatively the same background for NGAO as with LGSAO, 
even though the thermal background should be significantly lower for NGAO, due to the lower 
operating temperature of the NGAO system. We adopt the typical configuration used for stud¬ 
ies of gravitational lens systems (Fassnacht & Cohen 1998 [16]): K' filter, MCDS readout 
mode, read number A = 16, and narrow camera mode. A real observed PSF is used for the 
current AO system. The simulated PSF and performance characteristics of NGAO have been 
kindly provided by the NGAO team. 

4. The infrared imager and spectrograph IRIS (Larkin et al. 2010 [39]) working behind adaptive 
optics on the Thirty Meter Telescope (TMT) is selected to represent the performance of the 
extremely large telescopes that will be operational in the next decade. IRIS is expected to be 
close to diffraction limited in the K' band. The simulated PSF and performance characteristics 
of TMT-IRIS have been obtained from the IRIS team, from the TMT Exposure Time Calculator 
(ETC) ^ and from the paper by Do et al. (2014) [40]. 

5. Euclid is a space survey telescope planned to be launched by the European Space Agency. 
Euclid has two instruments: a high-resolution visible imager, and a NIR imaging spectrograph. 
We study the performance of the high-resolution visible imager, the more suitable to detailed 
gravitational lensing work, with the standard survey parameters. The instrument properties are 
taken from Schweitzer et al. 2010 [41], Penny et al. 2013 [42] and Cropper et al. 2014 [43]. 
Since the PSF is proprietary to the Euclid team, we adopt for simplicity a Gaussian function 
with full width half maximum O'.' 18. 

^Clearly, for survey telescopes such as LSST one can further increase the precision on y' by considering multiple bands 
at once, when they are available (Newton et al. 2011 [34], Auger et al. 2013 [35]). However, the multifilter analysis is 
beyond the scope of this paper, which is concerned with comparison across instruments/telescopes. 

^http://etc.stsci.edu 

^http://jwstetc.stsci.edu 

"^http: //tmt. mtk. nao. ac. j p/ETC_readme. html 
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6. We study the survey mode of the 2.4m WFIRST space mission, focusing on the 2000 square 
degree HLS Imaging Survey in the F184W filter. The instrument/survey properties are based on 
the WFIRST-Astrophysics Focused Telescope Assets Final Report (by the Science Definition 
Team and WFIRST Project, Spergel et al. 2013 [44]) and the WFIRST ETC. ^ Since the final 
PSF was not available at the time of this writing we adopt a Gaussian function with FWHM 
O'.'15. As an additional test, in order to investigate the sensitivity of our estimates to the choice 
of the PSF, we repeat the simulations assuming a Moffat (Moffat 1969 [45]) PSF, with different 
shape parameters. 

7. The Large Synoptic Survey Telescope (LSST) is designed to repeatedly image the entire visible 
Southern sky in six optical and near infrared bands (out to ~ 1 micron) every few nights for 
ten years. We take the survey and instrument properties from Ivezic et al. 2008 [46] and the 
LSST Science Requirements Document.^ The PSF of LSST will vary over time. We adopt a 
Gaussian PSF with FWHM O'.'? as representative of the stacked image quality of LSST in the 
i-band, and consider the full depth survey images in our analysis. 

We do not consider radio interferometers like ALMA or VLBI, given the radically different 
properties of the data, and the need for additional assumptions relating the optical/IR light to the 
radio emission. An assessment of their ability to provide sufficiently accurate high resolution images 
for cosmography is left for future work. 

3 The Lens Sample 

In this section we introduce our sample, taken to be illustrative of the kind of systems to be discovered 
in the next decade. 

3.1 Sample Selection 

We have chosen four prototypical systems for this exploration, characterized as faint or bright and 
double or quad depending on the photometry and image configuration. This set of four main choices 
covers regimes with different numbers of pixels above a certain S/N, which in turn depend on the 
number of images produced by the lens. 

The brighter lenses are selected to be representative of the majority of currently known lenses, 
selected from the Sloan Digital Sky Survey. For reference, the current largest quasar lens sample 
comes from the SDSS Quasar Lens Search (SQLS; Oguri et al. 2006 [47] ; Inada et al. 2012 [48]). The 
well-defined statistical quasar lens sample of it consists 26 lensed quasars which are selected from 
SDSS Data Release 7 (DR7) quasar catalog with Galactic extinction-corrected i-band magnitudes 
(15.0 < i < 19.1) in the low redshift (0.6 < z < 2.2), adopting morphological (image separations of 
1" < 0 < 20") and color (/-band magnitude difference between two images should be smaller than 
1.25 mag) selection algorithms. Among the 26 lensed quasars, the number of the quad and double 
images are 4 and 21 respectively, the remaining is a 5-image cluster lens. Meanwhile, 36 additional 
lensed quasars are identified from the DR7 quasar catalog, which contains 4 quad-image systems, 30 
double-image systems, the rest two are 5-image and 3-image cluster systems. 

The fainter lenses are selected to be representative of the fainter systems to be discovered in 
current and future surveys within the next decade (Oguri & Marshall 2010 [13]). We expect that the 
systems to be used to measure time delay distances in the next decade will span approximately the 

^http://wfirst-web.ipac.caltech.edu/w£Depc/w£Depc.j sp?etc 
^http: //WWW .Isst.org/files/docs/SRD.pdf 
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range covered by our sample. Note that some known lensed quasars are significantly brighter than 
any of the systems simulated here (e.g., Patnaik et al. 1992 [49]). Even though we expect more of 
those to be discovered in the near future and they will be useful, they will not represent the majority 
of the systems in future large samples, nor present any major challenge for follow-up, so we do not 
consider them here. 

To remain as realistic as possible, and take advantage of the high quality information available 
for the sample, the brighter mocks are based upon lenses in the Sloan Lens ACS Survey (SLACS; 
Bolton et al. 2006 [50], Treu et al. 2006 [51], Koopmans et al. 2006 [52], Gavazzi et al. 2007 [53], 
Bolton et al. 2008a [54], Gavazzi et al. 2008 [55], Bolton et al. 2008b [56], Treu et al. 2009 [57], 
Auger et al. 2009 [58], Auger et al. 2010 [59], Newton et al. 2011 [34], Shu et al. 2015 [60], Papers 
I-XII, respectively) and the fainter mocks on the Strong Lensing Legacy Survey (SL2S; More et al. 
2012 [61], Gavazzi et al. 2012 [62], Ruff et al. 2011 [63], Sonnenfeld et al. 2013 [64], Sonnenfeld 
et al. 2013 [65], Sonnenfeld et al. 2015 [66], Papers I-V, respectively, except for the first citation) 
samples. 

The sources are set at Zs = 1.071 and 2.77 and the deflectors are set at Zd = 0.351 and 0.783, 
respectively for the bright and faint systems. The redshifts are taken from the real systems used as an 
inspiration for this study, and are well within the range of the redshifts that we expect to find in future 
samples of lensed quasars (e.g. Oguri & Marshall 2010 [13]). A different choice within a realistic 
range would not have changed any of our conclusions. 

3.2 Parameters of the Sample of Mock Lenses 

The structural parameters of the four systems are listed in Table 2 and Table 3, which are used to 
describe the light model of the host and deflector galaxies and the mass model of the deflector. These 
models will be discussed in more detail in the next section. When not explicitly known from existing 
data, mock model parameters are assigned via plausibility arguments as specified below. The source 
and deflector of the bright lens system configurations are built along SLACS J0330-0020 as modelled 
by Bolton et al. 2008a (SLACS V, [54]), Auger et al. 2009 (SLACS IX, [58]), and Newton et al. 2011 
(SLACS XI, [34]). Lor t\\Q faint lens system configuration, the source and deflector have parameters 
from the Sonnenfeld et al. 2013 (SL2S III, [64]) model of SL2S J135949-r553550 (also, Sonnenfeld, 
A., 2015, private communication). 

Source and deflector magnitudes in the bright case for I and V bands are from Newton et al. 
2011 [34], for X-band are estimated based on the other colors, and typical spectral energy distribu¬ 
tions, i.e. an early-type galaxy for the deflector and a star forming galaxy for the source. Meanwhile, 
source magnitudes in ihQ faint case are arbitrarily set to 25.0 in all bands. Deflector magnitudes in 
Xht faint case are not directly available from the data, so they require some extrapolations. The value 
for /-band is based on Sonnenfeld et al. 2013 (SL2S III, [64]), while for K and V bands are estimated 
with K-corrections from typical spectral energy distributions. Unknown magnitudes in X-band are 
assigned from //-band ones via Kab - Hab- We note that the host galaxies of type-1 AGNs at these 
redshifts are generally brighter than the magnitudes adopted here (Bennert et al. 2011; [67]), so our 
assumptions are quite conservative. 

The effective (half-light) radii to be fitted in bright systems are evaluated separately in different 
bands, to get a precise estimate of the effective radius and its uncertainties. Differently, a unique 
effective radius is asserted across bandpasses for ihQ faint systems. Neglecting color gradients in this 
case does not affect the ability to recover the mass density profile in any significant way (Sonnenfeld 
et al. 2013 [64]). 

Lor both the faint and the bright systems, the source positions are assigned so as to map the 
source in either two or four images. This choice allows us to explore whether the number of images 
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makes any substantial difference while keeping the other properties of the systems fixed. 

In order to construct realistic systems for time delay measurements, point sources must be added 
to represent the lensed quasars. As described in the next section this is done in the image plane, in 
order to gain computational precision and efficiency. The magnification at the location of the quasar 
images is calculated by solving the lens equation using gravlens [68]. The source-plane magnitudes 
of the point sources are chosen to be somewhat brighter then the host galaxy, as it is typically the case 
for medium luminosity AGN (e.g., Bennert et al. 2011 [69]), and have realistic colors for AGN at 
the source redshits. In previous work (Peng et al. 2006a [70]; Peng et al. 2006b [71]; Rusu et al. 
2015 [33]), in general the lensed AGNs are found to be 0 to 3 magnitudes brighter than the host in 
the bands considered by our study, although there are a few exceptions where the host is brighter than 
the AGN. The adopted magnitudes are listed in Table 4. 

We note that the bright source configuration is too bright to be practical for TMT, given its 
sensitivity - the exposure time scales as where D is the telescope diameter, for background lim¬ 
ited point source exposures. Realistically, observations of such bright systems would be completely 
dominated by overheads related to target acquisition and will only be carried out in extremely rare 
circumstances. Therefore we do not simulate TMT observations for the bright systems. Even for 
the fainter system, we have chosen to make the AGN artificially fainter {K-26.()) in order to avoid 
saturating the central pixels while imaging deep enough to get a sufficient signal to noise ratio on the 
host galaxy. As will be shown below, the exposure time requirements for TMT are short even in these 
cases. 


4 Description of the inference process 

In this section, we summarize the process of modeling the mock lens systems. By comparing the 
distribution of the ratio between the inferred values and the input, we obtain the desired estimate of the 
statistical precision with which each parameter can be determined. As motivated in the introduction, 
we are interested primarily in the mass density profile slope y'. Therefore, we will only show the 
statistics for that parameter, even though all the parameters are varied simultaneously during the 
inference. We first describe the image generating process in Section 4.1, and then we describe the 
inference process in Section 4.2. 

4.1 Image generating process 

Before we fit the mock lens system we need to define the models to describe the source galaxy light, 
the lens galaxy light, and the lens mass. We start with the surface brightness distribution of the 
source and lens galaxies. The correspondence between image-plane and source-plane is given via the 
deflection angles, which in turn depend on the mass distribution of the lens. A PSF is used to convolve 
the light of lens galaxy, lensed source galaxy and point-source. Finally, we introduce counts noise, 
readout noise, and background noise. These steps are described in this subsection. 


4.1.1 Light Model 

We use the Sersic profile (Sersic 1963 [72], 1968 [73]) to describe the surface brightness profiles of 
the sources and the deflectors. The Sersic profile is described by 


1{R) = 4 exp 
R{x,y,q) ■■ 



^Jqx^ + y'^lq. 


(4.1) 

(4.2) 
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The amplitude Iq is the intensity at the effective radius the Sersic index n controls the degree of 
curvature of the radial light profile, the constant k is determined so that /^eff is the half light radius 
(Ciotti & Bertin 1999 [74]), and q denotes the axis ratio. Empirically, n increases with galaxy lumi¬ 
nosity, most galaxies being fitted by Sersic index in the range 0.6 < n < 10.0 (Merritt et al. 2006 
[75]). We adopt n-Afox the deflectors, typical of massive early-type galaxies, and n - 4/3 ovn- A 
for the sources, consistent with typical observed values for blue galaxies and AGN hosts. We checked 
by re-running simulations with different n for the source that our results do not depend significantly 
on this model assumption. 


4.1.2 Mass Model 

The lens mass profile is assigned within a class of power-law models. 


X(x,y) = Xc 


3 -/ 
1 + 




A-y 


Re 


(4.3) 


where is the axis ratio, y' is the radial power-law slope, the {X, Y} principal axes are rotated by 
the lens position angle w.r.to the canonical x (increasing to the West) and y (increasing to the North). 
Within this functional family, the deflections scale as R^~y with distance from the lens. The critical 
density is defined by 


c^Ds 

AjiGDaDas 


(4.4) 


in terms of the relative distances to the deflector (T>d), to the source (Ds), and between the deflector 
and the source (T>ds). The Einstein radius Re is chosen such that, in the spherical limit (q^ = 1), it 
encloses a mean surface density equal to This is also the radius of a ring traced by the host of the 
quasar when this is exactly aligned with the lens galaxy. 

Power-law models of elliptical galaxies (Evans 1994 [76]) have often been used with success 
over the years to model gravitational lenses. East methods to compute 2D deflections from elliptical 
power-law profiles have been given by Barkana 1998 [77], as a special case of the formalism by 
Schramm 1990 [78] for homoeoidal profiles. 

Lensed image profiles are obtained by inverse ray shooting. Eirst, each pixel position in the 
image plane is mapped back to the source plane via its corresponding deflection angle. Then, we 
exploit the fact that surface brightness is preserved by lensing to assign surface brightnesses from our 
putative source models to each detector pixel in the image plane. Einally, we convolve these images 
with the PSE appropriate for each instrument, yielding mock observations of the AGN host galaxy 
Einstein rings. 


4.1.3 Point Source Model 

In order to reduce computation time and complexity, the point source images are added to the light 
model by adding a PSE with appropriate normalization and position directly to the image plane. The 
normalization and position of the images for the center of the source in the source plane are computed 
for each configuration by solving the lens equation using gravlens [68]. 


4.1.4 Image Noise 

In the final step of our simulations we account for the effect of noise arising from the counts statistics, 
the background, and detector read out. As usual, the variance of noise per pixel is given by 


“ 'yjct -\- B ' t 


( 4 . 5 ) 








where C is the signal from clean lens system in electrons per second per pixel, t is the integration 
time in seconds, B is the sum of the sky background and detector dark current in electrons per second 
per pixel, A/read is the number of detector readouts, and R is the standard deviation of the read noise 
in electrons. Realistically, when long integration times are required, we set the number of readouts 
by requiring no single exposure be longer then a maximum time that depends on each instrument 
configuration. Considering the definite properties in different types of instruments, we pick 1000.0 
seconds as the crude maximum exposure time for space telescopes (HST and JWST), 300.0 sec¬ 
onds for the ground based telescopes (Keck and TMT) appropriately. Each of the surveys we model 
(Euclid, WEIRST, and LSST) have fixed maximum exposure times (590 seconds, 184 seconds, and 
15 seconds, respectively) which we use to compute the read noise per exposure. We then emulate 
the stacking of each survey’s exposures to reach the approximate expected total exposure time (2360, 
920 and 4500 seconds, for Euclid, WEIRST, and LSST, respectively). Examples of simulated images 
are shown in Eigure 3 to 6. 

4.2 Inferring the parameters 

The same code that is used to generate the mock Einstein Ring images is also used to fit them. The 
modelling procedure follows previous work done with real lens data (e.g.. Auger et al. 2011 [79], 
Auger et al. 2013 [35], Stark et al. 2013 [80]). We fit for the following non-linear parameters: posi¬ 
tion, half light radius, axis ratio of light, position angle, and Sersic index for the source and deflector 
brightness model (12 parameters); position, Einstein radius, the axis ratio of the mass distribution, 
position angle, and power-law slope for the deflector mass model (6 parameters); and QSO positions 
in the image plane (4 parameters for two-image lenses, 8 parameters for four-image systems). The 
sky level is fixed to the true value in this process. We note in particular that due to e.g., microlens- 
ing, millilensing, and differential reddening, the positions and fluxes of the QSO images are not tied 
to the lens model and they are therefore in effect treated as ‘foreground’ stars. Most priors on the 
non-linear parameters are taken to be uniform, but the prior on the position of the deflector mass is 
tied to the position of the deflector surface brightness; the mass is generally close to coincident with 
the light in galaxy-scale lenses. The fluxes/amplitudes of each of the surface brightness components 
(the source and deflector galaxies and the QSO images) are linear parameters and are determined by 
linear inversion for any given set of proposed non-linear paramerters. 

The goodness of fit is assessed through the standard image-plane;^^, by comparing the model 
and mock data surface brightness profiles. Uncertainties are obtained via MCMC exploration of 
the posterior, with the fluxes implicitly marginalized over. The likelihood is simply given by X oc 
exp[-^2/2]. 

The mock data images for different instruments are all 4" x 4", corresponding to the pixel 
numbers of 80 x 80 for HST, 130 x 130 for JWST, 400 x 400 for Keck (both LGSAO and NGAO), 
1000 X 1000 for TMT, 40 x 40 for Euclid, 36 x 36 for WEIRST, 20 x 20 for LSST approximately. The 
sampling numbers are fixed at 60000 for each inference progress. The acceptance ratios remain stable 
between 20% ~ 30%. Based on the parameters above mentioned, the inference is computationally 
expensive: typical run times range between 15 minutes and 4000 minutes per system on a linux 
desktop computer. In more specific terms, ~ 100 minutes for HST, ~ 150 minutes for JWST, ~ 300 
minutes for Keck (both LGSAO and NGAO), ~ 4000 minutes for TMT, ~ 40 minutes for Euclid and 
WEIRST, ~ 15 minutes for LSST. 

4.3 Exposure times 

The goal of this work is to investigate what combinations of instrument/telescope configuration and 
exposure time produce images of quality sufficient to determine y' with 0.02 precision. 
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For the telescopes/instruments operated in observatory mode, we are able to set the exposure 
time to any desired level. The minimum exposure time required is defined as the “target” expo¬ 
sure time. In addition, in order to quantify how the precision on 7 ' depends on exposure time, we 
also simulated images with 1/3 and 3x the target exposure time. This is meant to provide useful 
guidance for designing future experiments. However, for some combination of lens brightness and 
telescope/instrument configuration, the target exposure time is too short to be adopted in realistic 
observations. In practice, therefore these lenses will either not be observed with this setup or be ob¬ 
served for longer exposure time. Taking into account typical overheads for pointing and acquisition 
we evaluated that the brighter lenses are too bright to be observed with TMT. Likewise the brighter 
systems are likely too bright to be observed efficiently with JWST. We still simulated them, but set 
the minimum exposure time to 60s. Of course in some instances it may be beneficial to obtain longer 
exposures with JWST and TMT and thus exquisite data. 

For the survey telescopes, we simulated the planned exposure time as described above, and then 
simply considered the specific question of whether the ensuing data quality is sufficient to meet the 
requirement or not. 


5 Results 

Using the methods summarized in the previous section we generate 30 mock systems for most config¬ 
urations and 10 mock systems for TMT with the view of saving computing time. The mock systems 
have the same true parameters, but different noise realizations. Then we carry out the full inference 
for each one using the python Markov Chain Monte Carlo sampler pymc. The results are summarized 
in Table 5 and Figures 7 to 10. Table 6 summarizes the “target” exposure times for HST, JWST, 
Keck LGSAO and NGAO, TMT, i.e. those required in order to reach 0.02 uncertainty on 7 '. We 
reiterate that this is is only the random component of the uncertainty, and thus it represents a lower 
limit to the total uncertainty. For this reason we have set a rather stringent limit of 0.02, in order to 
leave room for additional uncertainties, such as, e.g., systematic errors related to PSF reconstruction 
and modeling errors. In practice, our target exposure times should be considered as minimum expo¬ 
sure times, i.e. necessary and not sufficient conditions. Assessing the contribution of those additional 
sources of uncertainty requires detailed modeling of specific systems and instruments (e.g., Suyu et 
al. 2014 [4]), and is left for future work. 

First of all, we find that HST can deliver the desired image quality with exposures of a few 
kiloseconds at most, consistent with published work (Suyu et al. 2013 [3]). Conversely, it seems that 
current AO capabilities - even setting aside the difficulties associated with the reconstruction of the 
PSF - are only sufficient to study the brighter systems, for reasonable exposure times. The improve¬ 
ment with NGAO will be substantial, cutting integration times down to HST-like for both the bright 
and the faint lenses. In this case, the higher background than from space is compensated by the ad¬ 
vantage of having smaller detector pixels, and thus better sampling of the PSF, and a larger telescope 
aperture. The NGAO performance may further improve if the thermal background is reduced with 
respect to LGSAO as planned and if the instrument has higher throughput than NIRC2. 

For TMT we concluded that 1200s of exposure will be more than sufficient for all systems 
considered here and therefore TMT is likely to be used either in “snapshot mode” or to take an image 
before taking a spectrum for determining the deflector velocity field and source redshift, or to follow¬ 
up fainter systems, yet to be discovered. Similary, JWST can obtain images of the desired quality in 
just a few minutes, and therefore it is likely to be used mostly to take images before taking a spectrum 
or to follow-up fainter systems. 
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For the surveys, the conclusion is that while the Euclid survey seems to have insufficient depth 
and resolution to meet our requirements on the targets considered here, LSST and WFIRST should 
be sufficient at least for the brighter ones. Of course, these surveys will cover vast fractions of the sky 
and will therefore provide useful images for all the brighter systems (Koopmans et al. 2009 [81]) at 
no additional cost. Those will be a great complement to the deeper pointed observation obtained with 
the other telescopes. WFIRST is planned to have a component of GO program, so it could meet the 
imaging requirements by integrating longer than in the survey mode. In fact, WFIRST is expected 
to have imaging performance superior to that of WFC3 on board HST, which has already been used 
to perform studies of time delay lenses, using integration times of a few orbits per system (~10ks; 
HST-GO-12889, PI: Suyu). Amongst the surveys considered here and modeled with a Gaussian PSF, 
WFIRST is the one with most potential for this application. Therefore, it is worth checking that our 
conclusion does not depend on the assumption of a Gaussian PSF In order to verify this, we have 
repeated our simulations adopting a more generic Moffat profile (Moffat 1969 [45]) 


yg-1 



(5.1) 


for the PSF, keeping the FWHM=2a V2F/^ - 1 fixed and varying the shape parameter y3. The results 
of the simulations are summarized in Table 7. The precision on y' decreases only slightly as the 
wings of PSF become more important (i.e. smaller j3), showing that our conclusions are robust to the 
choice of the PSF shape, as long as that is known. 


6 Summary 

Gravitational time delays are a powerful tool for cosmography. Transforming measured time delays 
into time distances requires modeling the gravitational potential of the main deflector with sufficient 
accuracy. Recent work has shown that this is possible provided that images of sufficient resolu¬ 
tion and signal to noise ratio are available. In anticipation of the hundreds to thousands of lensed 
quasars that are going to be discovered in the next decade from ground based low resolution imag¬ 
ing surveys (Oguri & Marshall 2010 [13]), we study the requirements for high resolution imaging 
follow-up. We consider a range of instrument/telescope configurations spanning from currently avail¬ 
able systems (HST/ACS and Keck/LGSAO/NIRC2) to planned observatories and surveys from space 
(JWST/NIRCAM, Euclid, WFIRST/HLS) and from the ground (Keck/NGAO/NIRC2, TMT/IRIS). 
In order to carry out a generic and systematic comparison across multiple telescopes we simulate four 
realistic lens systems, spanning the range of magnitudes and configuration expected for lenses to be 
discovered in the next decade. We also define a single metric, which is the ability to measure the 
slope of the mass density profile of an elliptical power law mass density profile. In order to guarantee 
that the images be sufficient to constrain the Fermat potential to a level so that the random uncer¬ 
tainties are subdominant with respect to other sources of uncertainty, we consider it necessary for 
them to contain enough information to constrain y' with 0.02 precision or better, which corresponds 
approximately to a 2% precision on time delay distance and hence Hq. We find that: 

• Our simulations show that HST can provide sufficient information with integration times of 
order 1-lOks, consistent with recent work based on real data. 

• Keck/LGSAO can provide sufficient information for the brighter systems to be discovered. The 
critical issue will be the accuracy with which the PSF can be reconstructed. Significant effort 
is under way to overcome this obstacle (Jolissaint et al. 2014 [82], Ragland et al. 2014 [83]). 
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• The planned Keck/NGAO system will improve the performance of Keck to make it comparable 
to that of HST, again provided that the PSF can be reconstructed with sufficient accuracy. 

• JWST/NIRCAM and TMT/IRIS will be able to deliver the desired image quality with short 
exposure times for any lens for which time delays are reasonably going to be available. Thus 
they will be ideal instruments for follow-up of time-delay lenses, especially if short images 
are followed by deeper spectroscopic observations to measure the redshift of the source and 
the spatially resolved velocity dispersion (and redshift) of the deflector. For TMT, and other 
ELTs, the abilty to reconstruct accurately the PSF of the AO system (Herriot et al. [84]) will 
be critical for time delay cosmography. 

• The survey modes of Euclid/WFIRST/LSST will likely be too shallow and/or have insufficient 
resolution except for the very brightest systems in the sky. Thus Euclid/WFIRST/LSST will 
make major contributions in the area of discoveries of these systems, but follow-up will be 
necessary for the fainter ones, for example with WFIRST in general observer mode. 

We conclude by emphasizing that our study is only concerned with precision, i.e. the study of 
random uncertainties. For this reason we have chosen a target precision of 0.02 on 7 ', that would 
be a subdominant contribution to the total error budget for the best currently studied systems. This 
is intended to leave room for additional sources of systematic uncertainty, like for example those 
related to the mass-sheet degeneracy. Within instrumental effects, an important term will stem from 
uncertainties in the knowledge of the PSF. Studies based on the comparison of multiple Hubble Space 
Telescope optical images have shown that in this case residual PSF uncertainties do not contribute 
more than the estimated total uncertainties. However, the contribution of the PSF to the total error 
budget may be more important for ground based Adaptive Optics data and future telescopes. Sig¬ 
nificant work is under way to characterize this source of uncertainty based on existing data (Rusu 
et al. 2015 [33]; Suyu et al. 2015, in preparation), even though we are not aware of a systematic 
investigation of the PSF-related uncertainties on the time-delay distance. We plan to carry out such a 
systematic investigation in the near future. 
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Table 1. Telescope and Instrument Properties 


Telescope 

Instrument 

Filter 

Zero Point 

Readout Noise 
(e"/pixel) 

Background Noise 
(e"/pixel/s) 

Pixel Size 
(arcsec) 

HST 

ACS 

E814W 

25.94 

4.20 

0.11 

0.050 

JWST 

NIRCAM 

E200W 

27.85 

9.00 

0.20 

0.032 

LGSAO 

NIRC2 

K' 

28.04 

5.75 

26 

0.010 

NGAO 

NIRC2 

K' 

28.04 

5.75 

26 

0.010 

TMT 

IRIS 

K' 

31.10 

2.00 

21 

0.004 

Euclid 

VIS 

r + i + z 

25.58 

4.50 

0.43 

0.100 

WEIRST 

— 

E184 

26.18 

5.00 

0.11 

0.110 

LSST 

— 

/ 

28.35 

5.00 

68 

0.200 


Observational facilities (telescopes and instruments) considered in this work. Zero points are given in the ABmag 
system. The VIS imager of Euclid spans the whole r + i + z wavelength range. The background and exposure time 
of Euclid will vary across the sky. Our fiducial values are representative of the average performance. The Keck 10-m 
telescope is considered both with the current (LGSAO) and next-generation (NGAO) adaptive optics capabilities. 


Table 2. Surface Brightness Profile Models 


Lens Name 

^efF 

(arcsec) 

q 

P.A. 

(deg) 

n 

Av 

(arcsec) 

Ay 

(arcsec) 

mi 

rriK 

mvis 

Mh 

Parameters for the source 

fainter system^ 

0.23 

0.92 

54.0 

4.0 

0.0662 

-0.167 

25.0 

25.0 

25.0 

25.0 

fainter system^ 

0.23 

0.92 

54.0 

4.0 

0.008 

0.298 

25.0 

25.0 

25.0 

25.0 

brighter system^ 

0.12 

0.77 

120.0 

1.33 

-0.195 

0.34 

22.73 

22.0 

23.46 

22.0 

brighter system^ 

0.12 

0.77 

120.0 

1.33 

0.01 

-0.005 

22.73 

22.0 

23.46 

22.0 

Parameters for the deflector 

fainter system 

1.76 

0.61 

-9.6 

4.0 

— 

— 

20.69 

19.7 

21.13 

19.7 

brighter system 

0.91 

0.81 

113.2 

4.0 

— 

— 

17.99 

16.5 

18.84 

16.5 


/?eff is the half light radius, q denotes the axis ratio. P.A. is with respect to the x-axis. The S^rsic index n controls 
the degree of curvature of the galaxy light profile. Magnitudes m are given in the ABmag system. The only difference 
among systems with sources the same brightness is in the source-position, which is set to produce either two or four 
images of the center. 

^ This configuration yields 4 QSO images. 

^ This configuration yields 2 QSO images. 


Table 3. Lens Model Parameters 


Lens Name 

z 

^Ein 

(arcsec) 

q 

P.A. 

(deg) 

r' 

fainter system 

0.783 

1.14 

0.6 

14.7 

2.0 

brighter system 

0.351 

1.1 

0.81 

113.2 

2.0 


Reiu is the Einstein radius, according to the definition by Kormann-h94. q denotes the axis ratio of the mass distribution. 
P.A. is anti-clockwise from the x-axis. The model deflector is a Singular Isothermal Ellipsoid (y=2), which is then fit 
using power-law models with variable y'. 
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Table 4. Point Source Parameters 



mag(I) 

mag(K) 

mag(VIS) 

mag(H) 

//(I) 

y^dl) 

//(III) 

//(IV) 

fainter system^ 

24.0 

24.0 

24.0 

24.0 

2.84 

3.75 

6.59 

2.78 

fainter system^ 

24.0 

24.0 

24.0 

24.0 

2.78 

2.48 

— 

— 

brighter system^ 

21.7 

21.5 

22.9 

21.5 

1.22 

6.10 

— 

— 

brighter system^ 

21.7 

21.5 

22.9 

21.5 

6.09 

8.22 

6.86 

7.53 


Magnitudes (in the ABmag system) and magnifications of multiple images. The source magnitude ‘mag’ is given 
in the source plane so that the observed magnitude in the image plane is given by mag - 2.5 log^oA^, where ji is the 
magnification of each image. We just consider TMT observations for the fainter system, and make the point source 
artificially fainter (^=26.0) in order to avoid saturating. 

^ This configuration yields 4 QSO images. 

^ This configuration yields 2 QSO images. 


Table 5. Estimated uncertainties on y' 




HST 

JWST 

LGSAO 

NGAO 

TMT 

Euclid 

WFIRST 

LSST 

fainter system^ 

l/3xtarget 

0.037 

0.048 

0.035 

0.033 

0.035 

— 

— 

— 


target 

0.019 

0.020 

0.022 

0.020 

0.021 

— 

— 

— 


3xtarget 

0.010 

0.010 

0.013 

0.011 

0.011 

— 

— 

— 


survey 

— 

— 

— — 

— 

0.034 

0.021 

0.038 

fainter system^ 

l/3xtarget 

0.034 

0.035 

0.042 

0.046 

0.034 

— 

— 

— 


target 

0.020 

0.019 

0.022 

0.022 

0.020 

— 

— 

— 


3xtarget 

0.011 

0.009 

0.014 

0.014 

0.009 

— 

— 

— 


survey 

— 

— 

— — 

— 

0.051 

0.023 

0.042 

brighter system^ 

l/3xtarget 

0.037 

0.013* 

0.052 

0.034 

— 

— 

— 

— 


target 

0.019 

0.006* 

0.022 

0.020 

— 

— 

— 

— 


3xtarget 

0.010 

0.004* 

0.009 

0.011 

— 

— 

— 

— 


survey 

— 

— 

— — 

— 

0.0107 

0.0044 

0.0025 

brighter system^ 

l/3xtarget 

0.042 

0.0087* 

0.042 

0.041 

— 

— 

— 

— 


target 

0.019 

0.0046* 

0.022 

0.020 

— 

— 

— 

— 


3xtarget 

0.011 

0.0032* 

0.011 

0.011 

— 

— 

— 

— 


survey 

— 

— 

— — 

— 

0.0036 

0.0016 

0.0007 


Estimated uncertainties on y' for each instrument, “target” represents target exposure time which is defined as that 
reaching 0.02 precision on the slope of the mass density profile of the mass model. * means the estimated uncertainties 
with artificial exposure time 60s, 180s, and 540s for JWST in brighter systems. 

^ This configuration yields 4 QSO images. 

^ This configuration yields 2 QSO images. 


Table 6. Exposure time requirements 


Instrument 

double 

faint 

bright 

quad 

faint 

bright 

HST 

6x lO^s 

360 s 

3 X lO^s 

150 s 

JWST 

690s 

180 s 

210s 

<60 s 

Keck (LGSAO) 

105 X lO^s 

3600 s 

75 X lO^s 

2400 s 

Keck (NGAO) 

18 X lO^s 

180 s 

12x lO^s 

150 s 

TMT 

1200s 

— 

1080s 

— 


Exposure times required to attain a 0.02 precision on the inferred y'. The Euclid, WFIRST and LSST surveys, are 
not quoted here because their exposure times are fixed. The default exposure times are sufficient to meet the 0.02 
requirement in the case of LSST and WFIRST, but not for Euclid. 





Table 7. Precision on y' with Gaussian and Moffat PSF for WFIRST 


Lens Name 

Gaussian 

Moffat 

08=1.5) 

Moffat 

08=3.0) 

II o 

Moffat 

08=6.0) 

fainter system^ 

0.021 

0.032 

0.028 

0.025 

0.023 

fainter system^ 

0.023 

0.036 

0.035 

0.033 

0.030 

brighter system^ 

0.0044 

0.0065 

0.0062 

0.0048 

0.0044 

brighter system^ 

0.0016 

0.0022 

0.0017 

0.0017 

0.0015 


Precision on the mass density profile slope y' with Gaussian and Moffat PSFs for WFIRST. This table represents the 
results for fainter and brighter systems with 4 and 2 QSO images. It is shown that as p increases the Moffat PSF tends 
to approximate the core of the Gaussian PSF. 

^ This configuration yields 4 QSO images. 

^ This configuration yields 2 QSO images. 
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Properties of simulated telescope/instrument configurations 
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Figure 1. Zero points in AB magnitudes for HST/ACS (blue), JWST/NIRCAM (green), Keck NIRC2 (assumed 
for both LGSAO and NGAO; red), TMT/IRIS (cyan), Euclid (magenta), WFIRST (yellow) and LSST (orange), 
corresponding to one count per second. The colored bars indicate the wavelength range of each instrumental 
setup considered in this work. 
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Figure 2. Montage of the point spread functions (PSFs) of each instrument. The upper row, from left to 
right, represents HST/ACS, Keck/NIRC2+LGSAO, Keck/NIRC2+NGAO, and JWST/NIRCAM. The lower 
row, from left to right, shows TMT/IRIS, Euclid, WFIRST, and LSST, respectively. Observed or simulated 
PSFs are used for HST, JWST, Keck (LGSAO & NGAO), and TMT. Gaussian PSFs are adopted for the other 
three survey instruments. 





Figure 3. Simulations of the fainter lens system with 4 QSO images. The simulated images are all 4" x 
4". The first 5 columns, from left to right, represent HST/ACS, Keck/NIRC2+LGSAO, Keck/NIRC2+NGAO, 
JWST/NIRCAM, and TMT; from top to bottom, the rows correspond to 1/3 x “target” exposure time, “target” 
exposure time, and 3 x “target” exposure time. “Target” exposure time is defined as the exposure time that 
yields 0.02 precision on the slope of the mass density profile of the mass model y'. The sixth column shows 
simulations of 3 surveys. From top to bottom we show LSST, Euclid, and WFIRST, with the default survey 
exposure times (4500s, 2360s, 920s, respectively). 


ass?’ 



















Figure 5. Simulations of the bright lens system with 4 QSO images. The simulated images are all 4" x 
4". The first 3 columns, from left to right, represent HST/ACS, Keck/NIRC2+LGSAO, Keck/NIRC2+NGAO; 
from top to bottom, the rows correspond to 1/3 x “target” exposure time, “target” exposure time, and 3 x 
“target” exposure time. “Target” exposure time is defined as the exposure time that yields 0.02 precision on 
the slope of the mass density profile of the mass model y'. The fourth column shows JWST simulations with 
3 fixed exposure times: 60, 180, 540 seconds, from top to bottom. The fifth column shows simulations of 3 
surveys. From top to bottom we show LSST, Euclid, and WFIRST, with the default survey exposure times 
(4500s, 2360s, 920s, respectively). TMT simulations are not shown since the system is considered too bright 
to be observed with this telescope in practice. 
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Figure 6. Same as Fig. 5, for the bright double imaged system. 
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Figure 7. Precision on the mass density profile slope y' as a function of exposure time. This figure shows the 
results for the fainter lens system with 4 QSO images. The input SIE mass slope is 2. The quantity is the 
output of the inference process. The histogram in each panel represents the distribution of inferred mass slope 
y' more directly. 
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Figure 8. Same as Fig. 7 for the fainter lens system with 2 images. 
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Figure 9. Same as Fig. 7 for the brighter lens system with 4 images. 


- 28 - 


















































Numbers o o m m m Numbers 


1.10 

1.05 


0.95 

0.90 


HST F814W 


■ t 

I . 

i 

4 

t = 120s 

t = 360s ■ 

t = 1080s 


1 ^ 



2.0 

2.5 3.0 

log t (s) 


3.5 



1.95 2.00 2.05 2.10 

Inferred Mass Slope 
LGSAO K' 



i t = 1200s 

1 J t = 3600s ■ 

1 ^ 4 t = 10800s 


I ♦ 


log t (s) 



1.10 

1.05 

I 1.00 

0.95 

0.90 



Inferred Mass Slope 
NGAO K' 





t = 60s 

- 


i 

t = 180s ■ 

1 I 


4 

t = 540s 

I i 

t 




log t (s) 


1.95 2.00 2.05 2.10 

Inferred Mass Slope 



1.85 1.90 1.95 2.00 2.05 2.10 

Inferred Mass Slope 


1.03 
1.02 
1.01 
I 1.00 

0.99 

0.98 

0.97 


Euclid VIS, WFIRST F184 & LSST I 


T 


i 

4 

Euclid: t = 2360s . 

WFIRST: t = 920s 

LSST: t = 4500s 

f ^ 


3.0 

3.5 

4.0 

4.5 


log t (s) 





1.85 1.90 1.95 2.00 2.05 2.10 

Inferred Mass Slope 


Figure 10. Same as Fig. 7 for the brighter lens system with 2 images. 
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